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ABSTRACT 


An approximate solution Is developed for the determination of tlie Inter- 
laminar normal and shear stresses In the vicinity of a crack In a three 
dimensional composite containing unidirectional linearly elastic fibers In 
an Infinite linearly elastic matrix. 

in order to reduce the complexity of the formulation, certain assump- 
tions are made as to the physically significant stresses to be retained. 

These simplifications reduce the partial differential equations of elasticity 
to differential-difference equations which are tractable using Fourier trans- 
form techniques. This "material modeling" approach Is in contrast with solu- 
tions developed by considering each lamina as a homogeneous, orthotiroplc layer. 
The resulting solution does not contain the classical singular stress field 
for the fibers and the influence of broken fibers on unbroken fibers is felt 
by a change in stress concentration factors. The matrix stresses however, 
are unbounded as the fiber spacing vanishes and an equivalent fiber-matrix 
geometry is proposed which gives the correct singular behavior. 

The numerical solution Is considered in detail and several specific ex- 
amples are presented. The potential for damaged or debonded zones to be gen- 
erated by an embedded crack is discussed, and stress concentration factors 
for fibers near the crack are given. Detailed comparisons are made between 
the present solution, the analogous two-dimensional problem, and corresponding 
shear- lag models. 
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INTRODUCTION 


In Investigations concerning the potential use of composite materials 
In advanced structures, conslderablif; attention has been given to problems 
associated with damage or flaws In laminates and to the resulting fatigue 
and fracture characteristics of the composite. The need for the capability 
to design damage tolerant composite structures with at least the same de- 
gree of confidence as now exists for metallic structures is essential If 
such materials are to be used In production. 

For composite laminates, the technique of hybrid material buffer strip 
embedded at regular Intervals throughout the laminate has been demonstrated 
for example see [ 1 1 , to have the capability of arresting through-the-lami - 
nate cracks, and in some cases to give a residual strength approaching the 
net section capacity of the laminate. The mechanism associated with the 
crack arrest and ultimate failure of such laminates has not, however, been 
explained nor have satisfactory methods been developed to allow for accu- 
rate prediction of crack arrest loads or ultimate failure in terms of the 
laminate geometry and materials. It does appear however, both f rom [ 1 ] and 
the experimental results reported In [2] that, for composite laminates con- 
taining embedded hybrid buffer strips or for single material laminates with 
weakened interface bonds, delamination between laminae appears to be direct 
iy associated with the laminates ability to restrict crack growth. 

The panel containing embedded buffer strips is considerably different 

than the related design where either bonded or riveted stringers are added 
externally without significantly interrupting the continuity of the basic 
panel. The panel and stringers then act much more independently during 
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crack growth and fracture than appears to be the case In composite laminates 
with Internal buffer strips. A detailed Investigation of the problem of 
riveted stringers reinforcing a crocked Isotropic elastic sheet Is presented 
In [ 3 ]. Similar studies for adhesively bonded external stringers are re- 
ported In [4] and [5]. An analogous study concerning composite laminates 
with internal buffer strips has not, to the knowledge of the authors, been 
completed. 

Considerable analytical work has been done for both Isotropic and ortho- 
tropic materials In Investigating the stresses In the vicinity of a crack 
tip for the crack being near, at, crossing, or along a material Interface. 
(See [ 6] , [ 7 ] , [ 8] , [ 9] , [ 10] , and [ 11].) These studies considered the 
laminate as a thin homogeneous Isotropic or orthotropic plane and Investi- 
gated the resulting two-dimensi'jnal stress state. Through the thickness 
variations In material properties and stresses were not considered. It 
appears that for the embedded hybrid buffer strip the significant stresses 
are, In fact, the interlaminar stresses and the three-dimensional problem 
must be considered. It does not seem reasonable to anticipate the develop- 
ment of a general solution in the manner of the above mentioned two-dimen- 
sional studies, as the extension to a three-dimensional strefs state with 
finite thickness would be exceedingly complicated. 

Some possibility of success might be had in modeling the lamlrate as 
in classical lamination theory, where, as discussed in [12] and [ I 3 ], the 
Interlaminar stresses near a free edge of a laminate have been approximated. 
Finite element or finite difference solutions also are a possibility; how- 
ever, without some Improvements in the current three-dimensional programs 
to account for through-the-thickness variations, the computation time seems 
prohibitbve. A different approach, wh I ch appears to have more potential, Is 



to model the laminate as a heterogeneous materia] consisting of separate 
fiber and matrix phases, each with Its own characteristic properties. If 
Che proper ossumptlof^', are made regarding the dominant stresses carried 
by each phase, the appropriate field equations can be simplified consider- 
ably while still retaining the fundamental nature of the solution. One 
major assumption usually made In such a model Is to neglect the normal 
stress In the matrix parallel to the fibers. The resulting solution then 
does not contain the classical singular stress field at the end of the 
notch but rather» the fibers feel the Influence of the size of the notch as 
a change In a stress concentration factor. 

Zweben In [ l4] presents a detailed discussion of this "material mod- 
eling" approach and In [ l4] and [ 15] a two-dimensional problem for a notched 
unidirectional laminate containing a damaged region with and without con- 
straint respectively, is considered. 

An earlier use of this technique was considered by Hedgepeth [ 16] and 
Hedgepeth and Van Dyke [17] for unidirectional fibers In a two-dimensional 
and three-dimensional array respectively. In each case however, the matrix 
was assumed to transmit only shear stresses and both studies reduced to 
one-d Imens Iona 1 solutions, Eringen and Kim, [ 18] extended the Hedgepeth 
model to Include normal stresses in the matrix and Investigated the two- 
dimensional problem of equally spaced parallel fibers In a single layer. 

The present paper extends the method [ 16], [ 17] , and [ 18] to the case 
of an arbitrary number of broken fibers in a three-dimensional composite 
containing equally spaced unidirectional fibers. The model developed below 
Is basically the same as [ l8] with some difference In the equilibrium equa- 
tion perpendicular to the fibers and the significant difference of the 
three-dimensional geometry. 
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Of course, the solution needed Is for o finite thickness laminate con- 
taining angle-ply laminae, each of possibly different material properties 
than the adjacent laminae as well as allowing for different materials with- 
in a particular lamina, It Is hoped that the present study will give some 
Insight into the behavior of Interlaminar stresses near a notch and aid In 
the development of more complete models. The extension of this Investiga- 
tion to Include angle-ply laminae Is currently being considered by the 
wrl ters. 
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FORMULA! I ON 


Consider an Infinite three-dimensional matrix containing parallel 
reinforcing fibers as Indlcar<jd In Figure (1), where both the fibers and 
the matrix are assumed to be linearly elastic. The development given 
below will be restricted to Identical parallel fibers with the spacing 
between fibers taken as equal. It Is assumed that the fibers carry 
only normal stresses and, In fact, support all the normal stress In the 
composite In the direction parallel to the fibers. The matrix supports 
shear stresses parallel to the fibers and normal stresses perpendicular 
to the fibers. Using the coordinate system given In Figure (1), a 
typical element can be Isolated as shown In Figure (2), where a particular 
fiber or matrix Interface Is Identified by the Indices m and n. Note 
that the composite Is doubly periodic In the y, z plane; however, this 
symmetry Is Initially not required of the solution and random fibers may 
be taken as broken with rhe resulting solution having no symmetry In 
the y, z plane. All fiber breaks are assumed, however, to occur on the 
plane x = 0 so symmetry does exist with respect to this plane. The 
reduction to the special case of broken fibers which are symmetric to 
the y, z axes is covered following the general development and a considerable 
reduction in computation time is realized. 

The equilibrium equations for the element shown in Figures (?,) and (3) are 


dff 


X 

m, n 


dx 


^ + h(T. 


xy 


m+l ,n 


xy 


} + K{t 


m,n 


xz 


■T. 


m,n+l 


xz 


} = 0 , 


( 1 ) 


|m,n 
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0 

^ m+ 1 , n 



m,n 


+ 


h 

2 


d r 

dx ‘V 


m+1 ,n 


+ 



0, and 


( 2 ) 


m,n+l 


m.n 


h d , 

2 dx ^^X 2 


m,n+I 


•t- X 


xz 


} = 0 . 


|m,n 


(3) 


The stress-strain and strain-displacements relations are approximated 
as follows, using .he geometry of Figures (A) and (5). 


du 

, = F ni.n 

xl dx 

lm,n 








(5) 




( 6 ) 


xy 


m+i ,n 


= xl ^ xl ■ ^7) 

M 2 dx m+1 ,n m,n h m+I ,n m,n 


T I = [w . , + w ] + ■^ [u . 1 " u ] } . 

xz M 2 dx r.i,n+l m,n h m,n+i m,n 

I m , n+ I 


(3) 


In the above equations u , v , and w are the x, y, and z components 

nifD rrijn m>n 

of the displacement of a point on the m,n Fiber located a distance x From the 

y,z plane. The stress a i Is the axial stress in the fiber and the 

^jm,n 


stresses a 


’ ^z\ ’ ^xv 

m,n ^im,n ^ 


m,n 


matrix. The constants Ep , Ej^^ and 


are normal and shear stresses in the 

m,n 

are Young's moduli of the fiber and the 
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matrix and the shear modulus of the matrix respectively, Ap Is the area 
of the* fiber, and h Is the shear transfer length. 

The boundary conditions are for an Infinite region, uniformly stressed 
at points remote to the mid-plane (y,z plane) containing an arbitrary number 
of broken fibers with the cracks always being at the mid-plane. Super- 
po^^'Icin may be used to consider two separate problems: the first, an 
Infinite region containing no broken fibers and having a uniform axial stress, 
0 ^, In the fibers, and the second, the infinite region with no stresses at 
Infinity and a compression stress of magnitude, , acting on the fibers 
at the cracks. The second problem will be considered below as It is the 
portion of the complete solution which Is non-trivlal. 

The boundary conditions are then 


u =v »w =0 forx=®, 

m,n m,n m,n 


(9) 


u =0 for all unbroken fiber at x = 0, (lO) 

m, n 


a 

x 


m,n 


-a 

o 


for all broken fibers at x = 0, 


(II) 



for al 1 f i bars at x = 0 


( 12 ) 



for all fibers atx=0 


( 13 ) 


The inclusion of the normal stresses in the matrix gives the necessary 
freedom to regu I re equations (.12) and (13). The above equations are 
normalized by letting 
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n “ x/h, U_ = u /h, V = V /h, W = w /h 
m,n m,n m,n m,n ’ m,n m,n 


(l^i) 


Y “ Ap/h^, Ej. = and Ej^ = E^/G 


M’ 


Using equations (M - (8) and equations {\k) , the equilibrium equations, 
equations (1), (2), and (3), may be written in terms of the displacements as 


2 dn ^Vl,n ’ '^m“l,n ^m,n+l " '^m,n-l^ 
^ ^m+l,n ^ Vl.n ■*■ ^m.n-H "-^m.n-l “ ''^m.n " ° ’ 


05) 


1 rl d^ 


+ V , ] + ^ {i- ^ [V + 2V + V , ] 

M m+l,n m,n m-l,n 2 2 dn^ m+I,n m,n m- 1 ,n 


+ ^ XI - U , ]} « 0 , and 

dn m+1,n m-l,n ’ 


(16) 


E,,[W ,”2W +W ("o' 'T~7 xi ■** 0 

M m,n+l m,n m,n"l 2 2 dq^ m,n+l m,n m,n-l 

+ [U - U ,]} = 0. 
dn m,n+l m,n-l 


07) 


It is now assumed that functions U (n, 9, <))), V (n, 9, *|>) , and 

W (hs 9. ft) exist such that the norma 1 i zed di spl acements U , V , and W 

m,n m,n m,n 

are the Fourier coefficients in a double Fourier series expansion. Then, 


CO 00 


U (n, 9, 4>) = s I u (n)e""'® e"'^‘^ 


r - ““ s = 


r,s 


08 ) 


CO CO 


V (n, 9, 4>) = ^ 


t V (n)e e and. 

r* 1 s 

|- S2 -00 5 =2 -no ' 


(19) 
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“(rS 


00 00 

W (n, 9, <&) = S E Wr 5 (n)e"'^'' e 

|« m mOO 5 » oOO * 


( 2 . 0 ) 


As the displacements are continuous functions of n, the representation given 
above Is necessarily valid. These equations may be Inverted to yield 



U (n, 9, (|i)e''’’® e*^*^ d 0 d(|), 


( 21 ) 


TT ir 

^ ^ d9d(Ii, and 


( 22 ) 


W (n, 9, <|i)e''"® d0d(J). 


(23) 


The equilibrium equations may then be written in terms of the functions 
U, V", and W as 


+ i sin (®) ^ ('i*^ ^ 

- 2 ] = 0 , 


(2i)) 


j [1 + cos (0)] 2E^ [1 - cos (0)] V + i sin (0) ^=0, and (25) 

y [1 + cos (<(>)] 0"- 2E^ [ I - cos (4.)] (7+ i sin (-(.) ^ = 0. (26) 



I ) 


These three second order differential equations can be written as uncoupled 
sixth order equations and If the boundary condition of vanishing displace- 
ments for large n» equation (9), is enforced, the solutions then will have 
the following form. 


U (n, 0, ‘t>) 
V (n, 0, <fi) 
W (pf 0, i|>) 


Aje + A 26 + Ase 

( 27 ) 

Bie'"’!" + + B36 "“s'' 

( 28 ) 

Cie’”!" + + Cje'^^n 

(29) 


where D^, D 2 , D 3 must have positive real parts. Substituting into equations 
( 25 ) and { 26 } the relationships between Cj^, and A|^» are determined as 


i sin (e) 

j [I + cos { 0 )] 0^ - 2¥j^ [1 - cos ( 9 )] 
i Aj^Dj^ sin ( 0 ) 

I" f' + (<|i)] 0^ - 2E^ [1 - cos 


( 30 ) 


(31) 


where k = I, 2, 3 . Equation (2^) then gives a sixth order algebraic equation 
for D|^ of the form 


aiD® + 02Dj^ + 0302 + ai4 = 0 . 


(32) 
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The coefficients Oj are 


yt p 

Ctl = -]J — [ 1 + cos (0)][ 1 + cos («)], 
a2“"2Y‘Ep'E^[ I - cos (0) cos 

“3 “ “h 0 1 - cos (0) ] {vEp£j_jr 1 “ cos + [ 1 + cos (tj)}] 

X [I - cos (e)]} + [ I - cos (<f>) ] (v'EpE'j^ [I - cos {0)] 

+ [ I + cos (0) ][ 1 - cos (i}i) ] , and 

at, “ 8^ [ I - cos (0)][ I “ cos (0)][cos (0) + cos (t()) - 2], 


Letting the six roots can be obtained and, in order to satisfy 

equation (9), the three roots with position real parts are taken. 

It Is noted that all roots of equation (32) vanish at 0 = ij) = 0 and 
unbounded roots exist at 0 “ ij) = +tt. The first is clearly seen from the 
fact that i f 0 = )(i = 0 then a2 “ 03 = at, = 0. For 0 = (|> = +ir, aj = a2 = 

0 and the resulting equation gives 

02 = + /S/yE- ' 0 = t|i <= IT. {37) 

1 1 


It also follows that the remaining roots must be unbounded and there- 
fore letting 0 = ({» = a ^ +ir, then 



This is the same value as the one singular root found in the two-dimensional 
case [ 18], An additional special case exist when either 9 or i|) equals zero 
as 04 is then zero and therefore D2 must vanish. The remaining two roots 
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are 


Dl, D3 


[ 2 E^[ l-cos(a)] 
[ l+cos(ct)] 




2 [ l+cos (g) ] 

-j 



( 39 ) 


where If 0 = 0» a * < 1 ) and for i{) = o» a » 6. The positive sign Is taken 
for Oi and the negative sign for D3. 

The functions A^, and A3, which depend on 0 and ^ , are next de- 
termined from the boundary conditions, where equation ( 9 ) Is satisfied by 
the proper choice of D^, Dg, D3. Equation (lO) gives 







t Ai + A2 + A3]e’’^®e*"‘*’d0d(j) = 0 


(i|0) 


for ail unbroken fibers, 


IT 


-It 


J 

-TT 


and equation (ii) gives 
[AiDi + A 0 D 2 + A3D3]e''^®e''^'*’d0dij) 


= 1 


( 1 * 1 ) 


for al I broken fibers . 

The conditions of vanishing shear stresses on n = 0 specified by 
equations (I2) and (13) require the following relations between the three 
functions A^ , A2, A3: 


A2 = -Av 




(i*2) 


{[ l+cos (iI))1D2^“ 1“cos (<ji) ] }{[ l+cos(e)]D 2 ^ " l-cos(0)]} 

{[ l+cos (<!>) ]Di^- 4 E|^[ l-cos((j))]}{[ 1+cos(0;]Di^ " 1|E|^[ 1 -cos (0) ] } 


and 


A3 = -Al 


Di^ - D2^ 


1^43) 


{[ l+cos ((J)) ] 03^ - 1*E^^[ 1-cos (([>) ]}{[ l+cos(e)]D3^ - l-cos(0)l> 


{[ l+cos (<j)) ]Di^ - ^E [ 1 -cos (<fi) ] }{[ l+cos(0)lDi- - 


l*Ej^[ 1 -cos (0) 1 } 


1^ 


Equations (^0) and {h\) then reduce to a set of dual inteqral equations In 
terms of the unknown function Proceeding In the same manner as 

[ l81, based on the orthogonality of the Fourier series, these dual Integral 
equations can be reduced to a set of simultaneous linear algebraic equations. 
Equation (^0) is Identically satisfied by assuming 


— [ Ai + A2 + A3] = Z Z 

r s 


-ire -isij) 


(^1*) 


where r and s correspond to the Indices of the broken fibers and the H 


r ,s 


are constants. Substituting from equations (^Z) and (^3) Into equations 
(40), (k\) , and then reduces equation (^1) to a set of algebraic equa- 

tions for the g. The set of equations are developed as follows. Refer- 
ring to equations i^2) and {^3), let 

A 2 ** 62^1 F 

A3 = R3A1 . (^6) 

Then, from equation (4^) 


— f 1 + 62 + PslAi = E E H e'"'® e“'®‘^ 

47t 2 r s 


m 


and equation (^’f) is then 
IT ir 


E E H. 
r s 


r ,s 


Di + S2D2 + 33D3 -ire -is^ ime Jn 4 . , 

e e e e d0d4i= 1 , 

J J ] + 02 ^3 

IT IT (^8) 


J 


where m, n, r, s correspond to broken fibers. With the constants H then 
known, the and Cj^ are specified and the stresses and displacement 

can be calculated. 

In the numerical evaluation of equation (it8 ) as well as the stress and 
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displacement equations, the singular behaviour of Di and O3 at 0 o ({, o + ir 
and the behaviour of Dj, D2, and D3 along the lines 0 or =» 0 or tt, must be 
accounted for. Noting the form of S2 ^nd 03 given by equations ( 1 | 0 ) , (lll), 
{h }) , and and the fact that Di= D3 ° 1+cos (a) ]' as 0 o ijj « a +n, 

It Is seen that 62 “ *^he limit. Proceeding to the limit as o -t- +7 t, 

the integrand of equation ( 48 ) then remains bounded and equals at 
e a (j) = TT, The analogous behavior occurs In the two-dimensional case [l8]. 
For either 0 or ^ equal to zero the Integrand of equation ( 48 ) Is bounded 
and vanishes at 0 = ij) » 0 . In a similar manner all stress and displacement 
equations have bounded Integrands at all points In the region of integration. 
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SYMMETRIC ARRAY OF BROKEN FIBERS 

The above equations can be simplified considerably for the special 
case of a symmetric array of broken fibers. For convenience in comparinq 
the two solutions only those equations which change due to the symmetry 
restriction will be recorded and they will be Indicated by the previous 
equation number with an asterisk. 

If a symmetric number of fibers is assumed to be broken with respect 

to the y and z axes It then follows that the displacement U (n) Is even 

m, n 

valued on the Indices m and n, while the displacement V (n) Is odd 

rn I n 

valued on m and even valued cn n, and the displacement W (n) is even 

m, n 

valued on m and odd valued on n. Equations (18) - (26) then become; 


u(n, 6, <l>) « El 

eo 

Si 

o 

It 

1_ 

s=o 

V(n, 0, <|)) = El 

O) 

Si 

r=l 

s=o 

w(n. 0i il*) = El 

CO 

Si 

o 

II 

S=1 


r,s 


r,s 


h ^ ^ - 

b „(n) = — f f U(n, 0, (}») cos(m6) cos(m))) d6d(t>, 
m « n 2 

tt^ 0 0 


_ 

V (n) = -T / / V(n, 0> 'i>) sin(m0) cos(n(j)) dedii, 

m.n 

’ 0 0 

4 ^ _ 

^ / W(n, 9, 4>) cos(m6) sin(n(f>) d 0 d 4 >, 

ni ♦ n *rr 

0 0 

Y^p + sin (9) ^ + sin((j)) ^ + 2V[ cos (9) + cos((l)}-2] 


(13)* 

(19) * 

( 20 ) * 

( 2 !)* 

( 22 )* 

(23)* 

0, (2^0* 



J7 


U 1 + cos(o)]^ 

A _t ^ 

dn*" 

i( I + cos4)liS 

^ dtl2 


- 2Ej[ I - cos(o)]V - sln(o) 
" ZT^ f ) -cos(<^)Jw - sfn(^) 


S!l 

dn 


a 0, 



(25)* 


( 26 )* 


Equations (27), (28), and (29) remain unchanged, and equations (30) 
and (31) become 




-A|^D|^ sln(e) 


“t I + cos(6)] - 2 eJ^ n “ cos(0)] 


and 


-A|^0|^ sln(<{i) 


M- cos(«<>)]d| “ ^ " c°s{^)] 


(30)* 


(30* 


The coefficients «i, «2» “3» ^^d are unchanged as are the roots to 

the characteristic equation ( 32 ). The remaining changes are In equations 
(^0), (M), m, m, and (1(8) as 

IT TI 

— / / [Ai + Ao + A3] cos(m0) cos (ntji) d0dc{i « 0, (l(0)* 

TT^ 0 0 

IT JT 

/ /[A^Dj + A2D2 + A3D3] cos(m 0 ) cos(n(|i) d0d<}i - I, ( 1(1 )* 

0 0 

-4-[Ai + A2 + A3] = E Z H cos(r 0 ) cos(siJ)), (A^f)* 

r 5 ’ ® 

if 

— [ ] + 62 + 63)''^! “2 Z H cos(r0) cos(si>), and (^7)* 

TT- r s 



)8 


r,s 


:t ’t 

/ / 


D1+62D2+83D3 

1+82+83 


cos(rO)cos(s(}l)cos(m0)cos(nt^) dOdf{» => 1. 


(itS)* 
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NUMERICAL SOLUTION AND RESULTS 


A computer code has been developed for the numerical solution of the 
above equations. The program was written In FORTRAN IV language for the 
NASA-Langley Research Center CDC computer.^ A maximum error of less than 
I percent was achieved with the average computation time being on the order 
of thirty minutes. 

The results presented below concern comparisons between the present 
three-dimensional solution, the two-dimensional studies of [16] and [18], 
and the three-dimensional shear-lag problem of [17]. Differences In the 
respective two- and three-dimensional analyses due to the Inclusion of the 
matrix norma) stresses In the model are central In the discussion. Detailed 
results are given apropos the manner of distribution of critical matrix nor- 
mal and shear stresses, and the rate of decay of fiber stress as a function 
of distance from broken fibers. All results are for “ 2.0, Fp = 5.2. 
Changes In material properties do not change the form of the solution and 
It was felt to be more Important to discuss the fundamental differences In 
the various solutions and to Investigate the significance of the geometric 
parameters, I.e. fiber spacing and number of broken fibers. 

It Is important to note that In (ll8], equation (2.1)), the equilibrium 
equation In the transverse direction, was written as 


0 

n 


+ 


b d / 
2 dy 


+ T 

n- 


) = 0 . 


^ The program is entitled, "Stresses in a Three-Dimensional Unidirectional 
Composite Containing Broken Fibers," and Is available through COSMIC, 

112 Barrow Hall, The University of Georgia, Athens, GA 30602, 
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This amounts to assuming that the fiber supports alt of the shear 
stress. We have changed this equation by replacing b with h to be consis- 
tent with the present analysis. The results for the two-dimensional geome- 
try presented In this paper were obtained by solving the appropriately 
modified equations of [ 18 ]. These changes gave Insignificant differences 
In the fiber stress but resulted In considerable changes In the matrix 
stresses. 

The constants y and h deserve further discussion as changes in para- 
meters can be Interpreted to Imply significantly different physical situa- 
tions. That Is, If ft Is taken as the distance between fiber boundaries, 
the relationship between y, h, and 6 and the aoproxlmatlons assumed for the 
matrix surface over which stresses are transferred and strains calculated 
are open to considerable freedom of definition, while still remaining within 
the validity of the formulation for both the t*wo- and three-dimensional 
models. The geometry of Figure 2 and the corresponding geometry of [l6], 

[17l» and [18] requires only that the distance h be consistent In all three 
equilibrium equations. 

For example, referring to Figures 2 and 3 and to equations (1) through 
(8), It is seen that h is the width in the vertical plane over which the 
she.,r stresses act in equation (l) , while In equations (2) and (3), it is 
the corresponding distance In the horizontal plane. In equations (5) through 
(8), it is related to the distance through which the matrix deformations take 
place, and should be bounded by the distance between fiber centers, b + 6, 
and the minimum spacing between fibers, 5. As seen below, the width in equa- 
tion {!) need not be the same as the distance in equations (2) and (3) or (5) 
through (8) . 


In Figures 6 and 7 three different geometries and corresponding assump- 
tions Yi h, and 5 arc given for the two- and three-dimensional models, respec- 
tively, The different value of y foi' each case Is due to the particular geome- 
try assumed and Is derived by writing the equilibrium equation In the axial 
direction, accounting for the specific fiber shape and shear transfer region. 
Table 1 gives values of y for particular values of 6 for each case of Figures 
6 and 7, 

Case I in Figures 6 and 7 corresponds to assuming that the fibers have 
the same material properties In shear and transverse normal directions as the 
matrix, and would be more appropriate for nearly equal fiber-matrix properties. 

Cases 11 and III are more realistic for a typical composite, In which tne 
fibers usually have a much higher modulus than the matrix, and make the assump- 
tion that the normal and shear strains occur over a critical distance equal to 
the minimum distance between fibers. The main difference between II and Ml Is 
the assumed shape of the fiber. The decision as to whether a square or circular 
fiber cross-section is more appropriate is based on considering the solutions 
for close fiber spacing and comparing the manner in which the stresses increase 
with decreasing ^ with previous elasticity solutions [19] and [20], In these 
studies the stresses between two circular cylinders under the action of inplane 
loading and transverse shear, respectively, were considered. In both studies 
the numerical results Indicated a l/v'5” behavior for the maximum stresses between 
rigid cylinders as the distance .S approached zero. For elastic cylinders un- 
bounded stresses did not occur. Cases II and III are equivalent to assuming 
infinite transverse normal and shear material properties for the fibers as all 
the strain Is assumed to occur over the distance between fibers; therefore. It 
seems consistent to seek a model having matrix stresses which behave as i/vT" 
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if possible. From Figures 6 and 7» Cases II and III, the matrix material 
supporting the stresses decreases like 1/h and l/h^, respectively, and the 
stresses should Increase more repldly with decreasing h In III than II. 

By considering the present solution along with [|6], [17] » find [18] and 
expandlig the equations for small values of h, where h = 6, the matrix 
stresses for all Four models are found to Increase as 0(l//h~") and 0(]/h) 
for Cases II and III, respectively. 

Case II Is then more appropriate for close fiber spacing as It agrees 
with the continuum solutions [I 9 ] and [20], and it is suggested that the 
actual fiber cross-section be represented as an "equivalent-' square cross- 
section having equal areas. All results given below are for Case II with 
the corresponding y, h, and 6 relationship of Figures 6 and 7 and Table 1. 

Figure 8 gives the value of the maximum stresses as a function of h for 
a square array of twenty five broken fibers, (five broken fibers In the two- 
dimensional study [18]). The stresses for [I 8 ] are larger for all values of 
spacing h due to the added constraint of the present geometry. As the number 
of broken fibers increase, the two solutions should approach each other. 

This is indicated in Table 2 for the particular case of one by eleven broken 
fibers i.e. (0,~5) through (0,5) where the stresses adjacent to the (l,0) 
fiber are compared with the stresses adjacent to the first unbroken fiber for 
one broken fiber In [18]. 

Both the normal and shear stresses are of considerable importance in 
formulating a failure criterion capable of predicting matrix damage. It is 
seen in Figures 9 and 10 that compressive normal stresses exist over part of 
the matrix region between the last broken and the first unbroken fiber, both 
in [ 18 ] and the present study. The shear stresses decrease due to the differ- 
ent boundary conditions at n = 0, from the corresponding shear-lag stresses; 
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therefore, it seems likely that the shear-lag model would lead to a prediction 
of greater matrix damage parallel to the fibers than the present model. In- 
deed, related two-dimensional shear-lag studies concerning matrix splitting 
for a unidirectional laminate (see [21] for one broken fiber and [22] for an 
arbitrary number of broken fibers with splitting and yielding) indicate that 
only a very small increase In the stress required to initiate damage can be 
tolerated before unstable splitting occurs. The splitting develops in the 
region adjacent to the last broken fiber which, as mentioned above, also has 
compressive stresses using the model of [18] without damage. The analysis of 
[22] is being extended to include a model having norma! stresses in the 
matrix, and It Is hoped that the results will be more consistent with experi- 
mental evidence, tn which case local, stable, splitting has been observed [23] 
and [2^]. Such behavior has also been observed at the NASA-Langley Research 
Center for Boron/Aluminum composites by C. C, Poe, (private communication). 

An equally interesting observation is noted for the region between the 
first and second unbroken fibers. As indicated in Figures 9 and 10, the shear 
stress is considerably reduced from its value on the opposite side of the 
fiber; however, it is larger than the corresponding shear-lag solution. More 
importantly, It is opposite sign than the shear-lag stress. The normal stress 
is tensile and more than fifty percent larger than the maximum shear stress 
in the adjacent region. First-ply matrix splitting has been noted for bonded 
joints [25] which are very similar geometrically to the region in the vicinity 
of the notch in the present work and the large tensile stress found here cer- 
tainly admits the possibility of such damage. 

Table 3 gives results using the present solution for a square array of 
broken fibers as well as analogous results from [l6], [17], and [l8]. The 
maximum fiber stress occurs on the plane of the break (n = 0) and in the first 
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unbroken fiber at the center of the sides for the three-dimensional problems 
and in the first unbroken fiber in the two-dimensional problems. For both 
the two- and three-dimensional solutions, the maximum fiber stress is seen 
to be slightly smaller for the shear-lag model. Both the fiber and matrix 
stresses are larger in the two-dimensional case [18] than the present solution. 

Figure II gives the magnitude of the maximum fiber stress in the first 
unbroken fiber as a function of the number of broken fibers for a square 
array, and compares these values with [16], [17], and [18]. Also, for the 
one particular case of eighty-one (nine) broken fibers, the stress in the 
first four unbroken fibers Is given. The decrease in stress is largest be- 
tween the first and second intact fibers with the rate of change being rela- 
tively small from that point on. As seen in the figure, the influence of the 
broken fibers decreases more slowly for the two-dimensional model [l8] than 
for the present three-dimensional model. 



CONCLUSIONS 


The significant difference between the present model (either two- or 
three-dimensional) and the corresponding shear-lag model Is in the change 
in the manner of distribution of the critical shear stress in the matrix 
and the introduction of matrix normal stresses. As noted, the critical 
shear stresses are reduced on the order of fifty percent over the correspon* 
ding shear-lag solution and the matrix norma] stresses are found to be of 
sufficient magnitude to suggest potential matrix damage. For close fiber 
spacing both the normal and shear stresses in the matrix are shown to in- 
crease without bound and an equivalent cross-section is proposed which gives 
the correct singular behavior based on previously published elasticity solu- 
tions . 

Surprisingly, with these major changes in matrix stresses, the fiber 
stresses are relatively unchanged, i.e., approximately three percent higher 
than the corresponding shear-lag solution. The fiber stress in all solutions 
is essentially independent of fiber spacing. 

It then seems reasonable to expect, that any attempt to develop an 
analysis for the Investigation of matrix damage must be based on a more com- 
plete model than the shear-lag assumption. 
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^Two-dimensional geometry, Figure 6, 
^Three-dimensional geometry, Figure 7* 

TABLE t. Relationship between the fiber spacing <S and the 
constant y for the geometry of Figures 6 and 7* 
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Model 

(3-0)1 

(2-D)^ 

(3-0)3 

r — 

(2-D)*^ 

0 

X 

1.371 

1.372 

1.327 

1 .333 

oy 

0.316 

0.316 

- 

- 

’xy 

-0.238 

- 0.238 

-0.3l<3 

- 0 . 3 ^^ 


^ present solution, fiber O,0). 

^ two-dimensional solution [I8], fiber (1). 

^ three-dimensional shear-lag [17]» fiber (1,0). 
two-dimensional shear-lag [16], fiber (I). 


TABLE 2. Comparison of maximum stresses for an array of 

one by eleven broken fibers in the three-dimensional 
model and one broken fiber in the two-dimensional 
case. The geometry of Figures 6 and 7 case II with 
h=6=1.0, Y=l,0ts used. 
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M» 


a (2-0)3 

X 

j (3-0)“ 

(2-0)“ 

1 

1.171 

1.38k 

I.lk6 

1.333 

9 

1.521 

1.908 

l.k56 

1 .828 

25 

1.821 

2.318 

1.728 

2.216 

'♦S 

2.087 

2.666 

1.967 

2.5k6 

81 

2.321 

2.97k 

2.181 

2.838 


^ M equals the total number of broken fibers In the three- 
dimensional case. for the two-dimensional case) 

^ present solution. 

^ two-dimensional solution [18]. 

** shear-lag solutions [16], (17). 


TABLE 3. Maximum fiber stress (stress concentration factor) for 
a square array of broken fibers, h ■ 5 • 2.0, y ■ 0.5. 
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FIGURE I. Three-dimensional array of parallel fibers. 



Figure 2. Free-body diagram for a typical element of the (m,n) fiber 
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Plane "n" where (a, b, c) and (a*, b', c') are the points before 
and after 'j^format ion. j ^ 

Note: To Insure sytwnetry Y,« t Iv , + v ] rather than 

^ ‘ 2 dx m+l ,n n,n 

T“ V 
dx m,n 


Figure *4. Displacements for a region containing the m,n fiber. 
(x,y plane) 



Plane "n" «here (b, c. d) and (b', c", d") are three points before 
and after deformation. I ^ 

Mote: To insure svmmetrv — [w w ] rather than 

^ Z dx m, n+ I m, n 

dx 


Figure 5. Oisplacements for a region containing the m,n fiber. 
(x ,2 plane) 
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Fiber stress for a square array of broken fibers compared with [l 6 ], ll 7 l, and 1 ) 8 ] 


